function plotFit(min_x, max_x, mu, sigma, theta, p)
    hold on;
    x = (min_x - 15:0.05:max_x + 25)';

    X_poly = polyFeatures(x, p);
    size(X_poly)
    X_poly = (X_poly - repmat(mu, size(X_poly, 1), 1)) ./ repmat(sigma, size(X_poly, 1), 1);
    size(X_poly)

    X_poly = [ones(size(X_poly, 1), 1), X_poly];
    size(X_poly)
    plot(x, X_poly * theta, '--', 'LineWidth', 2);
    hold off;

end
